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Abstract 

Quantum canonical transformations are used to derive the integral repre- 
sentations and Kummer solutions of the confluent hyper geometric and hyperge- 
ometric equations. Integral representations of the solutions of the non-periodic 
three body Toda equation are also found. The derivation of these representa- 
tions motivate the form of a two-dimensional generalized hypergeometric equa- 
tion which contains the non-periodic Toda equation as a special case and whose 
solutions may be obtained by quantum canonical transformation. 

1 Introduction 

The confluent hypergeometric and hypergeometric equations underlie many of the 
special functions of mathematical physics and, by extension, many of the exactly 
soluble problems of quantum mechanics in one dimension. The solutions of these 
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equations have of course been well understood for more than a hundred years (see 
e.g. This paper returns to derive the integral representations and Kummer so- 
lutions of the confluent hypergeometric and hypergeometric equations using a new 
approach whose motivation comes from physics. The purpose of this is two-fold: to 
give an elegant and economical derivation of time-honored results and to illustrate 
the use of some basic tools which have application to solving more difficult problems. 
As an illustration, simple derivations of two integral representations of the wavefunc- 
tions for the (non-periodic) three-body Toda equation are given. This motivates a 
2-dimensional generalized hypergeometric equation which includes the non-periodic 
Toda equation as a special case. 

The origin of the method to be used lies in quantum mechanics. Given a linear 
differential equation, it can be interpreted as a time-independent Schrodinger equa- 
tion at a given energy, where the dependent variable is taken to be the wavefunction 
and the differential operator to be the Hamiltonian in the coordinate representation. 
(For convenience, the discussion will be given assuming the differential equation is in 
terms of a single variable; the generalization to many variables is straightforward.) 
The differential operator cum Hamiltonian can be abstracted from the coordinate 
representation and viewed directly as a function of non-commuting variables x and 
d which satisfy the commutation relation [x,d] = —1. The coordinate representation 
consists in taking d = d/dx. We use d instead of the physicists' p to avoid unnecessary 
factors of i; expressions in terms of p are simply obtained by substituting d = ip. 

The method consists in making a quantum canonical transformation from the 
given Hamiltonian to another one having known solutions. A quantum canonical 
transformation is a change of variables (x, d) i— > (x'(x,d),d'(x,d)) which preserves 
the algebraic relation [x,d] = — 1 = [x',d'\. It is so called by analogy to a classi- 
cal canonical transformation which is a change of the phase space variables which 
preserves the Poisson bracket. Quantum canonical transformation are produced by 
acting with a function C(x, d), 

q^q' = CqC-\ d ^ & = CdC~ 1 . (1) 

Under this transformation, any function of x and d transforms 

H(x, d) i-> H'{x, d) = CH(x, d)C~ l = H{CxC~\ CdC~ l ). (2) 

A solution ip of Hip = is obtained from a solution ip' of CHC^ip' = H'lp' = by 

i> = C- 1 ^'. (3) 
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Care must be taken if the kernel of C as an operator in the coordinate representation 
is non-trivial as not all ip obtained in this way will be solutions of H (see [||). 

2 Quantum canonical transformations 

To be more precise about the quantum canonical transformations, one must address 
the nature of functions C(x, d) and their representation as operators. By introducing 

and d^ 1 as the algebraic inverses of x and d (xx~ l = x~ x x = 1, dd^ 1 = d~ 1 d = 1), 
one can define an algebra U consisting of functions of the variables x and d and 
their inverses, consistent with the commutation relation [x, d] = 1. The algebra 
is constructed so that every element C € U has an algebraic inverse C _1 G U, 
and so that elements transform under a quantum canonical transformation by (0). 
Details are given in Ref. ||. The precise class of functions has not been identified 
with a named class, but it essentially consists of those functions and their algebraic 
inverses which do not involve distributions in themselves or any of their derivatives 
(the functions may have poles, branch points or essential singularities). Every linear 
differential operator with suitably smooth coefficients corresponds to an element of 
this algebra. Very likely it is possible to generalize the analysis to treat operators 
whose coefficients are not smooth, but this has not been investigated. 

The functions C(x, d) and their inverses are well-defined as elements of the algebra 
IA. Their realization as operators in the coordinate representation is more subtle. (For 
convenience, the same notation will be used for the elements of the algebra and their 
realization as operators.) Care must be taken because the operators corresponding 
to C and C^ 1 are not strictly inverses because their kernels may be non-vanishing. 
One can define them as inverses when restricted to act on suitable subspaces with 
the kernels projected out@, but it is practical to think of them as being inverses up 
to arbitrary linear combinations of elements in their respective kernels. For example, 
the inverse of differentiation, d = d/dx, is indefinite integration, d^ 1 = j x dx, up to 
an arbitrary additive constant (an element of ker d) . 

One might be concerned about possible problems in defining the operator re- 
alization of the quantum canonical transformations. The property that makes the 
quantum canonical transformations useful is that they can be decomposed into a 
product of elementary transformations, each of which has a well-defined action on 
wavefunctions|2|]. There are three elementary canonical transformations: 
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1) similarity transformations, C = exp(J x f(q)dq), 

/X 
f(q)dq)ip(x), (4) 

2) point canonical transformations, Pf( x ), 

x^f(x), d ^J7^ d: p f(x)^(x) =^{f{x)), (5) 
and 3) interchange, /, 

1 r°° 

x^-id, d^-ix: Iip(x) = ^ 1/2 J_^dye lxy tp(y). (6) 

The functions in each transformation are allowed to be complex. In the differential 
equation context, these correspond to similarity transformations, change of variables, 
and Fourier transform. By composing them, one obtains ordered functions of opera- 
tors whose behavior on wavefunctions is well-defined. In particular, 

/ix pQ 
f(q)dq)r 1 = exp( J f(q)dq) 

and 

IPf{ix)I~ l = Pf(d) 

give the analogous similarity (x, d) i— > (x + f(d),d) and point canonical transforma- 
tions (x,d) i— > (f'(d)~ 1 x,f(d)) involving d. These are the basic tools for making 
quantum canonical transformations. They encode in an efficient way the details of 
performing the standard manipulations of differential equations. 

As an illustration, consider the transformation generated by t xd which will be used 
later. This can be decomposed as 

t xd = exp(lntx<9) = P lnx exp(lntd)P e *. (7) 

The middle transformation is recognized as a translation while the outer ones are 
point canonical transformations. Their collective action on the algebraic element x 
is then 

t xd xt -xd = p lnxexp (\ nt d)p exX p lnxexp (-\ nt d)p eX ( 8 ) 
= Pin x exp(x + In t)P e x 
= tx, 
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and their action on a wavefunction is also to rescale x 

t x9 ijj(x) = P lnx exp{]ntd)P e ,'ip(x) = ip{tx). (9) 
Alternatively, one could simply have observed that since 

g, — l^xdg, ^x~ 1 xdx ^xd+l 

one has 

t xd xt -xd = xt xd+i t -xd = xt 

For simple manipulations like this, it is often not necessary to explicitly decompose 
the transformation into elementary canonical transformations, though it may be done. 

There is one main difference between the use of quantum canonical transforma- 
tions for physics and for simply solving differential equations. In physics, one often 
wants to find a transformation which takes the entire spectrum of an operator into 
that of another. In particular, one would like to find a canonical transformation from 
an interacting to a free theory. If one is interested instead in an equation like the 
hypergeometric equation, one is not concerned with the spectrum of the operator, but 
rather with a particular solution. The significance of this is that one does not need to 
transform to an operator for which one can find the general solution, only to one for 
which a single solution is easily identified. This will generate a single solution of the 
operator one is interested in. Other solutions can be obtained from other canonical 
transformations, especially from those which act as symmetries of one's operator, 
taking it into another operator of the same form but with different parameters. 



3 Confluent Hypergeometric Function 

The differential operator for the confluent hypergeometric equation can be repre- 
sented by the Hamiltonian 

H = xd 2 + (b - x)d - a. (10) 

The confluent hypergeometric functions are the solutions ip of 

Hip = 0. (11) 

The Hamiltonian is easily canonically transformed to an equation soluble on inspec- 
tion after it is rewritten as 

H = (b + xd)d - (a + xd). (12) 
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The quantum canonical transformation produced by 

T(a + xo) 

transforms the derivative 

d » CdC- 1 = p^Q. (14) 
b + xd 

This is easily verified by observing 

dT(a + xd)d~ 1 = T(a + dx) = T(a + xd + 1). 

Since C commutes with b + xd and a + xd, it transforms the Hamiltonian to 

H ia) = CHC' 1 = {a + xd)(d- 1). (15) 

One of the solutions of f/w^w = is clearly ipw = e x . The complementary solution 
is easily found but won't be needed. (In a problem for which one was interested in 
the complete spectrum Hip = Eip, one wouldn't be satisfied with this form. One 
would look for further transformations leading to, say, a free Hamiltonian H' = d 2 .) 
The solution of the confluent hypergeometric equation which corresponds to ip( a ' = 

€> X IS 

* = CTV W = S^x^e"- (16) 
F(b + xo) 

Expanding the exponential in a power series in x and allowing the operators to act, 
one finds the familiar power series representation of the confluent hypergeometric 
function 

^ T(a + n) x n 

* = (17) 

= ^M(a,b,x). 

One can obtain an integral representation of fllTD by recognizing that C^ 1 times 
T(b — a) is a beta function, 

r(o+ ,f )r( a 6 - o) = fdu^ii -tr«- 1 (is) 

T{b + xd) Jo K J v ' 
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> 9?(a) > 0). Applying this to e x , and using (|9|), one finds 

-M(a,b,x)= dtt a -\l-t) h - a - l e tx . (19) 



r(6) 
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This is a familiar integral representation from which others may be obtained by 
change of variables. 

A Barnes-type contour integral representation is also easily obtained by expressing 

e x as 

1 MOO 

e* = — / dsT(-s)(-xy. (20) 
Applying C~ l to this, one has 

r ( a ) w/ j \ 1 r , Vt-s)V(a + s), xs 

i m w i »=»L 4 W" i) ' (21) 

One must choose the contour of integration so that it separates the poles of r(— s) 
and r(a + s). Again, this is the standard result [||]. 

An expression for the other familiar confluent hypergeometric functional U (a, b, x) 
can be found from a second canonical transformation. The transformation 

° = T(a + xd)T(l-b-xd) ^ 

transforms the momentum 

8«-2±f9fl. (23) 
b + xd 

The transformed Hamiltonian is 

#("') = (a + xd)(-d- 1). (24) 

This has the immediate solution e~ x . The confluent hypergeometric function corre- 
sponding to e~ x is 

I»r(l + a - b)U(a, 6, x) = T(a + xd)T(l -b- xd)e~ x . (25) 

One recognizes another instance of the beta function, 

dtt a+xd ^{l + tf- a -\ (26) 



r(a + xd)T(l-b-xd) _ r, + + a +xd -ln , 



r(l + a - b) 
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Applying this toe 1 gives 

POO 

T(a)U(a,b,x) = dtt^Hl + tf-^e^K (27) 
Jo 

Additional forms of solutions can be found by applying symmetry transformations 
of the equation. These give the Kummer relations between solutions. A symmetry is 
a canonical transformation which takes an operator into another operator of the same 
form but with different parameters. Since we are interested here in equations of the 
form Hip = 0, the term "symmetry" will be broadened to include transformations 
which produce an operator of the form of H up to an overall factor, which can 
be divided out. Let H(a,b) denote the confluent hypergeometric operator ([10|). A 
symmetry C gives 

CH(a,b)C- 1 = FH(a',b') (28) 

where F is an overall factor. Given a solution ip' of H(a', b'),ip = C~ Xr i\)' is a solution of 
H(a, b) (assuming here that the kernels of F and C vanish; a more general discussion 
is possible^). 

The simplest symmetries are those obtained by similarity transformation. Since 
the coefficient of the highest derivative term is x, the transformation C\ = x~ a , taking 

d^ddCV 1 = d+-, (29) 

x 

will produce a linear derivative term of correct form. The transformed Hamiltonian 
is 

ff(«) = x d 2 + (2a + b-x)d+ a - a{l - ^ - (a + a). (30) 

x 

Taking a = or 1 — b cancels the undesired 1/x potential term, and leaves an operator 
of confluent hypergeometric form. The symmetry is 

dHfafyCi 1 = H(a + a,b + 2a), (31) 

where a = 0, 1 — b. The transformation with a = is trivial, but the other produces 
a solution of H(a, b) 

^ = CfV (a) = x l ~ b M(a - b + 1, 2 - b,x). (32) 

Since this solution has different behavior than (|17|) at the regular singular point x — 0, 
it is an independent solution. 
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A second symmetry transformation is obtained from the transformation C% = 

e -/3x 

d^d + (3. (33) 

This leads to 

ffW = xd 2 + {b-{l- 2(3)x)d + b(3-a + ((3 2 - /3)x. (34) 

Requiring the coefficient of x to vanish gives (3=1. The transformed Hamiltonian is 

#W = xd 2 + (b + x)d + b-a. (35) 

This is not of confluent hypergeometric form because the sign of x in the linear term 
is wrong. A second transformation S = exp(i7rxd) 

x 1 — > —x, d 1 — > —9 (36) 

takes ii/"^ to in canonical form (PH|) up to an overall sign. This sign is not 
important because the equation to be satisfied is H(a, b)ip = and an overall factor 
can be divided out. The symmetry here is 

SC 2 H(a, b)C 2 l S- 1 = -H(b -a,b). (37) 

A solution of flTUD is then given by 

ijj = C^S- 1 ^ = e x M(b - a, b, -x). (38) 

Combining these two transformations gives a third transformation, 

C t SC 2 H{a, b)C 2 l S- l C^ = -H(b -a + a,b + 2a), (39) 

where a = 1 — b. This leads to the solution 

^j = x 1 - b e x M(l-a,2-b, -x). (40) 

Inspection of the behavior near the regular singular point x = shows that 

M(a,b,x) = e x M(b- a,b,-x), (41) 
x 1 - b M(a-b + 1,2 -b,x) = x 1_6 e !e M(l - a, 2 - b, -x). 

These are the Kummer transformations 0. The above symmetry transformations 
may also be applied to U(a, b, x) to find similar relations. 
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4 Hypergeometric Equation 

A similar story can be repeated for the hypergeometric equation. The Hamiltonian 
corresponding to the hypergeometric equation is[|J 

H = x{l- x)d 2 + (c-(a + b+ l)x)d - ab. (42) 

This can be rewritten as 

H = (c + xd)d - (a + xd)(b + xd). (43) 
The canonical transformation 

C = r(c + ^ (44) 

T(a + xd)T(b + xd) V ; 

transforms the derivative 

d CdC^ = ^^){b^xd) 

c + xd K ' 

The transformed Hamiltonian is 

H {a) = (a + xd)(b + xd)(d -I). (46) 

This has ip^ = e x as a solution. The solution of H corresponding to ip^ is the 
familiar hypergeometric function 

^HpF(a,6;c;x) = C"V (47) 

T(a + xd)T(b + xd) 



e x 



T(c + xd) 
T(a + n)T(b + n) x n 
F(c + n) n\ 



E 



Using (pCf ), one obtains a Mellin-Barnes integral representation |4]] 

«F(, , c; ,) = -L r r(a + y + .)r(-.) 

1 (CI Z7TZ J-ioo 1 (C + SI 
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where the contour separates the poles of T(a + s)T(b + s) from those of r(— s), a, b, c 
are not negative integers, and a — b is not an integer. 

Another integral representation is found by recognizing that C _1 times 
r(c — b) contains a beta function. Thus, one finds 

T(b)T(c-b) _ T(b + xd)T(c-b)T(a + xd) 

r(c) na,*,c,x) - T{c + xd) r(a) e W 

j i ib-\-xd — 1 / 1 .\c-b-l ^( a + X &) x 
{I t) 

(9?(c) > > 0). This can be further simplified by expanding e x in a power series 
and allowing T(a + x<9) to act on it. This gives 

r(a) ^ r(a)r(n + i) x • lbUj 

Using the gamma function inversion formula, r(z)T(l — z) = 7rcsc7r,2, this gives 
na + xd) _ f r(l-q) 

= (1-*)-. 

Inserting this in ( ^9j) , and letting t xd act on (1 — x)~ a , one obtains the familiar integral 
represent at ion 

r(6)r(c-6) F = f\ n b-U 1 _ t y-b-u 1 _ tx) -a (52) 

r(c) jo 

The set of Kummer's 24 solutions of the hypergeometric equation can be obtained 
by symmetry transformations which take the hypergeometric equation into itself. Let 
H(a,b,c) represent the hypergeometric operator fl4"2|) . Since a quadratic function of 
the coordinate multiplies the d 2 term, the most general shift of the derivative which 
produces a linear term of the correct form is generated by C\ = x~ a (l — x) -13 , 

d h-> ddcr 1 = d+-- -J—. (53) 
x 1 — x 

This produces a new Hamiltonian 

if to = x (l - x)d 2 + (c + 2a - (a + b + 1 + 2a + 2d)x)d (54) 

a 2 - (1 - c)a f3 2 + (a + b- c)8 , _ w _ 

+ — + \ — - (a + (3 + a) (a + (3 + b). 

x 1 — x 
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This is in hypergeometric form if the undesired potential terms are cancelled by taking 

a 2 - (1 - c)a = (55) 
(3 2 + (a + b - c)(3 = 0. 

The solutions of the conditions are a = 0, 1 — c and (3 = 0, c — a — b, and these give 
three transformed hypergeometric functions. (Note that these are the exponents of 
the regular singular points at and 1.) The symmetry transformation is 

CxH(a, b, c)Cf 1 = H(a + a + (3, b + a + (3, c + 2a). (56) 

A solution of the original hypergeometric equation is given in terms of the transformed 
one as 

$ = CrV (a) = x a (l - xfF{a + (3 + a,a + p + b;c + 2a;x) (57) 
For the four choices of (a,j3), one finds 

(58) 





= F{a,b; 


c; 


x) 




= x^Fi 


a 


— c + 1, b — c + 1; 2 — c; x) 




= (1-x 


r 


" a - b F(c-a,c-b; c;x) 




= x 1 -^! 




x y- a - b F(l -b,l-a;2-c;x 



Inspecting the behavior at the regular singular point x — 0, one finds that = ip^ 
and t/>( 2 ) = ip^ . For 1, if)^ and tp^ are independent solutions. 

In addition to similarity transformations, one can make point canonical transfor- 
mations. The transformation Py x takes 

x i— ► 1/x, d t— > — x 2 d. (59) 

This transforms the Hamiltonian to 

#(°) = - x [x(l - x)d 2 + (l-a-6-(2- c)x)d + — ]. (60) 

This is not in canonical form because of the term in 1/x. This can be cancelled by 
making the similarity transformation C\ fl53|). The transformed Hamiltonian is 



= -x[x(l - x)d 2 + (1 - a - b + 2a - (2a + 2(3 + 2 - c)x)d (61) 
+ a 2 - a(a + b) + ab {3 2 + /3(a + b - c) 



x 1 — x 

-{a + (3){l - c + a + (3)}. 
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Taking a = a,b and (3 = 0, c — a — b cancels the undesired potentials and returns 
this to the hypergeometric form. (Note that a is the exponent of the regular singular 
point at oo.) One has 

C x P ljx E{a, b, c)P 1/x C{ 1 = -xH(a + (3,1 - c + a + (3,1 - a - b + 2a). (62) 

The overall factor of x is irrelevant since the equation is H^ip® = 0, and it may be 
divided out. The solutions of the original hypergeometric are equation 

jj, = P 1/x Cf V b) = x~ a - p {x - lfF{a + (3, 1 - c + a + (3; 1 - a - b + 2a; l/x) (63) 

For the different choices of {a, (3), one finds 

^,(5) = x -*F(a,l-c + a;l + a-b;l/x), (64) 

^,(6) = x -bF(b,l-c + b;l+b-a;l/x), 

^) = x b - c (x-l) c - a - b F(c-b,l-b;l+a-b;l/x), 

^ (8) = x a - c {x-l) c - a - b F{c-a,l-a;l+b-a;l/x). 

Checking behavior at the regular singular point x = oo, one finds = and 

^(6) 

A second point canonical transformation P\- x also leads to new solutions. This 
transforms 

x i — 1 — x, d i — ► —d. (65) 

The transformed Hamiltonian is 

#(«) = x (l - x)d 2 + (l + a + b-c-(a + b + l)x)d - ab. (66) 

This is in hypergeometric form. Applying the transformation C\, with a = 0,c — a — b, 
(3 = 0, 1 — c (because c in the original hypergeometric equation has been replaced by 
1 + a + b — c) , one obtains 

ifW = x (i_ a; ) ( 9 2 + (i + a + 6- c + 2a-(2a + 2/3 + a + 6 + l)a;)9 

-(a + (3 + a)(a + (3 + b). (67) 

One has 

C 1 P 1 _ x H(a, b, c]P x _ x C^ x = H(a + a + (3,b + a + (3,1 + a + b - c + 2a), (68) 
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with a = 0, c — a — b, (3 = 0, 1 — c. The solutions to the original hypergeometric 
equation are 

ip = Pi_ x C7fV (6) =x f3 (l- x) a F(a + p + a,a + p + b;l + a + b-c + 2a;l-x). (69) 
For the different choices of (a, P), these are 

(70) 



^(9) 


= F{ 


a,b; 


1 + a + b — c; 


1 - x), 






^(10) 


= (1 


— x] 


f- a - b F(c - b, 


c — a;l + c — a — 


M 


-x), 




= x 1 


~ C F( 


a — c + 1, b — 


c+ l;l + a + 6- 


c;l 


- ar), 




= x 1 




- x) c - a - b F(l 


— 6, 1 — a; 1 + c — 


- a — 


b;l — x 



Checking behavior at the regular singular point x — 1, one finds t^ 9 - 1 = ?/>( n ) and 

^(10) = ^(12)_ 

The previous two transformations can be composed to lead to further solutions. 
First, one can apply C\P\/ X C\P\- X . One must remember that C\ depends on two 
parameters a, (3 that are not reflected in the notation. These parameters are deter- 
mined by the form of the hypergeometric function that C\ acts upon, and are in 
general different for each factor of C\. Using (|S^) and (|62|) , one finds 

C x P xjx C x P^ x H{a,b,c)P^ x C^P xjx C^ = (71) 
= CiP 1/x H(a + a 1 + pi,b + a 1 + Pi,l + a + b-c + 2a 1 )P 1/x C x 1 
= —xHia.2 + @2, c — a — b — 2a x + a2 + P2, 1 — a — b — 2a i — 2(3\ + 2a 2 )) 

where 

a x = 0, c — a — b, (72) 

Px = 0, 1-c, 

«2 = ctx + Pi + a, a x + Px + b, 

p 2 = 0, 1-c- 2ft. 

There are only four distinct forms of solution, and without loss of generality one can 
take (ax,Px) = (0,0). This gives 

V> (13) = (1 -x)- a F(a,c-b-l + a-b;—^—), (73) 

1 — x 

V> (14) = (1 -x)- b F{b,c-a-l + b-a; ' 



1 — x' 
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^ (15) = {-l)- a x 1 - c (x - l) c - a - x F(l + a-c,l -b;l + a-b;-^— ), 

1 — x 

^ (16) = (-l)-V~ c (a;-l) c - 6 - 1 F(l + 6-c,l-a;l + 6-a;-^). 

X 37 

A second transformation is induced by C2 = C\P\- x C\P\j x . Without loss of 
generality, one can take («2,ft) = (0,0), this makes the symmetry transformation 

C 2 H(a, b, c)C^ = (x - l)H(an + ft, 1 - c + a ± + ft, 1 + a + b - c + 2ft), (74) 

where a\ = a, b, ft = 0, c — a — b. The solutions of the hypergeometric equation are 

^( 17 ) = x - a F(a,a-c+l;a + b-c+l;l--), (75) 

37 

^,(18) = x -6F(6,6-c+l;a + 6-c+l;l - -), 

V> (19) = x b - c (a: - l) c - a - fe F(c - 6, 1 — 6; c - a — 6 + 1; 1 ), 

^(20) = x «-c( a: _ 1 )c-o-6 F ( c _ 0jl _ a;c _ a _ 6+1;1 __) j 

Finally, the transformation induced by 

C3 — CiPi- x C\P\/ x C-]_P\- X i 
taking (cn^ft) = (0,0) and (a 2 ,ft) = (a, 0), makes gives the symmetry 

C 3 H(a, b, c)Cz l = {x- l)H(a + a 3 + ft, c - b + a 3 + ft, c + 2a 3 ), (76) 
where 0:3 = 0, 1 — c, ft = 0, b — a. This gives the hypergeometric solutions 

^ (21) = (l-x)- a F(a,c-b;c;^-), (77) 

^(22) = (_ 1 )l-c x l-c( 1 _ x )c-a-l jF ( a _ c+1>1 _ 6;2 _ c; _^.) j 

^ (23) = (l-x)- fe F(6,c-a;c;^-), 

rr — 1 

^,(24) = (_ 1 )l-c x l-c( 1 _ x )c-6-l F ( 6 _ c+1>1 _ a;2 _ c; _^.)_ 

37 X 

This completes the construction of the 24 Kummer solutions of the hypergeometric 
equation. 
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5 Three-body Toda Quantum Mechanics 



As an application of canonical transformations to a more difficult problem, consider 
the nonperiodic three-body Toda potential. The Schrodinger operator is 

H = -~(d? + d\ + dl) + e qi ~ q2 + e 92 " 93 . (78) 

This operator does not separate in any of the standard orthogonal coordinate sys- 
tems, yet integral representations for its eigenfunctions can be constructed. The 
n-body Toda potential, generalizing ( [T§D in the obvious way, is physically interesting 
as an exactly solvable nonlinear oscillator which may be obtained from the Korteweg- 
deVries equation by an appropriate discretization 0. 

Two integral representations of the three-body Toda wavefunctions are quoted in 
the review by Olshanetsky and Perelomov || , and both can be economically derived 
using canonical transformations. The result of Vinogradov and Takhtadjan was 
originally found using methods from number theory. A second form was obtained 
in || by constructing power series solutions for the wavefunctions and then by sim- 
ply stating the associated integral representations. The derivations of these integral 
representations by canonical transformation are analogous to the solutions of the 
hypergeometric equations, in that a sequence of transformations, one involving the 
Gamma function, is used to transform the equation to a form where a solution is obvi- 
ous by inspection. An alternative solution in which the Toda problem is transformed 
to a free one with vanishing potential is discussed in || . 

The result of Vinogradov and Takhtadjan J/J for the 3-body Toda wavefunction is 
the most easily obtained and so we shall begin with it. The 3-body Toda Hamiltonian 
(ff8|) can be transformed into center-of-mass coordinates with separated interaction 
potential by the transformation P c .o.m. 

91 i-> §9i + §92 + 1<?3, <9i i \di + d 3 

92 >-> -§9i + hh + ks, d 2 ^ -ld 1 + \d 2 + dz (79) 

93 >-> -f 9i - |?2 + p3, <9 3 ~\d 2 + d 3 . 

The inverse transformation M is 

(91,92,93) >-» (§9i - §?2,|?2 - §93,9i + 92 + 9s) (80) 

(9 1; 5 3 ) ^ (|^ - §<9 2 - §9 3 , 1^ + p 2 - p 3 , \d x + \d 2 + 19 3 ). 

The resulting Hamiltonian takes the form 

H a = - l -(dl + d\ - dxdz) - hi + e 291 + e 2 " 2 . (81) 
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The following sequence of transformations then reduces the Hamiltonian to a form 
for which an eigenfunction can be found by inspection. 

[VT1] exp(-ln2(<9i + c> 2 )) : 
[VT2] exp{-ik( gi - q 2 )) : 
[VT3] r ^ d 1 +3ik ^ d 2 -Sik ^ 



qi h-> 


9i - In 2, 


9 a H 


- & 




qi i-> 


<?2 - In 2, 


<9 2 H 


-> <9 2 




q\ i-> 






- 3l 


+ i& 


<?2 l-> 




<9 2 h 


-> <9 2 


— zA; 


e 2 9l ,_ 


9i+3ife„2<7i 




- <9i 




e 2q 2 ,_ 


&2-3ik 2q? 
£>i+<9 2 e ' 


<9 2 i- 


-> <9 2 





(82) 

Suppressing the free-particle kinetic term — |6f describing the motion of the center 
of mass, the Hamiltonian transforms as follows: 



H a = -\(dl + d\ - d x d 2 ) + e 2(?1 + e 2 " 2 

[VT1] i-> -\{dl + dl-d x d 2 -e 2 <K -e 2 ^) 

[VT2] i-> -i(d 2 + 3ikd 1 + d 2 - 3ikd 2 - d 1 d 2 - e 2qi - e 2q2 ) + ^ 

[VT3] - --^ [{dl + 3 zk)(d 2 -e^) + (d 2 -3zk)(d 2 -e^)] + ^ 



33) 



Inspection shows that this final Hamiltonian has B v {e qx )C v {e q2 ) as an eigenfunction 
where B u , C u stand for modified Bessel functions, each either I u , K u or a linear 
combination of them. The eigenvalue of this eigenfunction is \{—v 2 + 3k 2 ). 

The eigenfunction of the Toda Hamiltonian ( ff8| ) (with zero center-of-mass mo- 
mentum) corresponding to this eigenfunction is 

Vv(?i>92) = N k,vPc.l. M . ex P ((#i + 30 ln2)exp (ik(qi - g 2 )) (84) 

p/ di +3ik \T^{ d2—3ik \ 



2 



Recognizing the product of Gamma functions as a Beta function, and using a familiar 
integral representation for the Beta function, one finds the result 



= ^,^c"oMe ifc(9l - 92) rdtt- 1 -^ (92 - 3ife) (l + t)i (ai+02) ^(2e f ' 1 )C(2e^) 

1 poo 3 

= N k , v e2 ik ^- 2q2+q ^ / dtt~ 1+ 2 ik (85) 



&,(2*/(l + t)e2^-^))C(2 v /(l + ^e^ 2 "^) 
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The physically interesting solution is the wavefunction which dies in the classically 
forbidden region, and this is found by taking B v and C v to be K„. With the identifi- 
cations k = — \\{t — s), v = — 1 + | (3s + 1), this agrees with the result of Vinogradov 
and Takhtadjan, as quoted in ||. 

The solution given by Bruschi et al. || can similarly be found by this approach. 
The key idea is to begin by factoring out the asymptotic plane-wave behaviour. De- 
note the momentum of this asymptotic plane wave by A = (Ai,A 2 ,A 3 ) in the co- 
ordinate system where the Toda Hamiltonian takes the form (|78|). Also let a = 
Ai — A2, b = A2 — A3. The sequence of transformations to a theory for which an 
eigenfunction can be recognized by inspection is 

H= -l(dl + dl + dl) + ew + ew 
e- iXq : » -i((9 1 +iA 1 ) 2 + (9 2 + iA 2 ) 2 + (9 3 + iA 3 ) 2 ) 

P c o m : 1 * - \ {d\ + d\ - d 1 d 2 + 2iad 1 + 2ibd 2 - 4e 2 ^ - 4e 2< ? 2 ) 
-f d 2 3 - i(Ai + A 2 + A 3 )«9 3 + \{\\ + A 2 + A 2 ) 

where P c .o.m. as in (|79|). Suppressing the eigenvalue ^(Af +A|+A§) and terms involving 
<9 3 which describe the center of mass motion, the transformations continue as 

e -ln2(d 1+ d 2 ) . ^ _|( 9 2 + Q2 _ di Q 2 + 2iaQi + 2{hd2 _ e 2 Ql _ e 2q 2 j 
T(l+ia+|ai) 



r(-i6-|a 2 )r(i+i(a+6)+f9i+|a 2 ) 



^(^felW^ 1 + 2ia )^i + 2i ^ a + b ^ ~ e 91 ) (87) 
+ (<9 2 + 2i&)(<9 2 + 2i(a + 6)<9 2 + e 2 « 2 )] 

e i(a+6)( gi + 92 ) . _> __L_[( 9i + _ 6) )( 9 2 _ e 2 9l + (a + 6)2) 

+(<9 2 + i(6 - a)){dj + e 2 <^ + (a + 6) 2 )] 

The final Hamiltonian in this sequence has i^a+f,) (e gi ) H^ +b Je 92 ) as an eigenfunc- 
tion with eigenvalue |(A 2 + A 2 , + A3). This eigenfunction will lead to a solution which 
is asymptotically free in the classically allowed region with asymptotic momentum 
A = (Ai,A 2 ,A 3 ) and which dies in the classically forbidden region. Other solutions 
to the equation can be obtained by using other products of Bessel functions. This 
corresponds to the Toda solution of Bruschi et al. Assembling the full transformation, 
this solution is 

*b(a) - n ^p-i > 2 f^2) r ( ~ ib ~ ^ 2 )r(i + i(a + b) + |gi + |gp 
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The product of Gamma functions is once again a Beta function and has the integral 
representation (with parameters at the limits of convergence) 

r(-i6-i^)r(i + i(o + 6) + i ^+i^) = 



r(i + ia + 



= e- b ^ r d~z 3 (i - Ssy^-^C'^ 92 - 



The translation operator exp(i7r<9 2 /2) acts to transform H$ +b Je q2 ) into K^ a+b )(e q2 ) 
times some constant factors. 

Letting the exp[ln2(<9i + d 2 )] factor act, one finds 

Ml) = — ^Pol m r ^ 3 (1 - 5 3 )- 1 - ia -^ 1 C 1_i6 "^ (90) 

7T JO 

e -i(a + 6)( 9l+ , 2 ) e (a + 6)^. (a+b)(2e9l)K . (a+b)(2e92)> 

The integral representation for the modified Bessel function 

r°° / e 2q \ 

J dxx"- 1 exp [ifi(x - —)) = 2e uq e lun/2 K_ l/ (2fie q ), (91) 

(Im(fj,) > 0, Im(e 2q fi) < 0) for v = — i(a + 6), with /x = 1 and e 9 appropriately 
displaced into the complex plane for convergence, can be used to replace the Bessel 
functions. Acting with the translation operators appearing in the Beta function, the 
result is 



1 V. poo 

2rr Jo 

/ dz 1 / dz 2 (2 1 ^ 2 )- i ( a+ ^- 1 e i(zi+Z2) exp ( - i— - i ). 

JO JO V Z 2 {1 - Z 3 ) 7 



Finally, absorbing the constant factors into the normalization and making the change 
of variables z 3 = z 3 /(ziz 2 ), one reaches the form of Bruschi et al. 



00 00 00 



MQ) = N' x e lX - q l / / d Zl dz 2 dz 3 ( Zl z 2 - Z3 )K^i)-i4 a 3-a 2 )-i (Q3) 

,i(*i+* a ) exp C _ i £l e ii'Q2 _ i^gm-my 

ZiZ 2 - Z 3 ^3 
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Note that this corrects some typos in the formula quoted in U, most importantly 
concerning the range of integration. Note also that this is an eigenfunction of (|78|) 
and not of the Hamiltonian of the form discussed in Section 12 of |J. 

The power of the quantum canonical transformations is evident in the simplicity 
of these derivations. 



6 2-d Generalized Hypergeometric Functions 

The example of the 3-body Toda problem suggests a nontrivial generalization of the 
hypergeometric equation to two dimensions. The center of mass form of the Toda 



problem (81) has the form of two Bessel operators in q\ and q2 with a cross term 
8182- The Gamma function transformation [VT3] in fl82|) removes the cross term 
and effectively separates the problem so that a product of two Bessel functions is a 
solution. In the resulting full solution, the Gamma function transformation acts as 
a convolution of the Bessel functions. Since a Bessel function is essentially a special 
form of hypergeometric function, it is natural to consider the function obtained by 
similarly convolving two (generalized) hypergeometric functions. 

The characteristic feature of generalized hypergeometric functions is the Gamma 
function structure of their series expansions 

,F.K... >ar ;fc fc;,) = E |^^ g ("> 

^0 I15=i r (% ) 111=1 r iPk + n) n\ 



-F s is a solution of the operator 



H=[]l(b k +xd))d-]l(a j +xd), (95) 
fc=i 3=1 



The canonical transformation 



transforms H to 



r _ nu nh + xd) 



H a =[H(a j +xd))(d-l), (97) 



20 



which has e x solution. 

In a two-dimensional form, the intent is to couple two hypergeometric equations 
with a cross term. Define 

A x = flia^+x^), (98) 
ji=i 

A 2 = f[(a%+x 2 d 2 ). 

32=1 

A natural equation is 

H = A x dx + A 2 d 2 - [(x^i) 2 + {x 2 d 2 f - x 1 d 1 x 2 d 2 (99) 
+a(x 1 d 1 - x 2 d 2 ) + 0\. 

By only making transformations which depend on X\di and x 2 d 2 , only the isolated 
di and d 2 factors will transform. The transformation 

C 1 = T(x 1 d 1 +x 2 d 2 ) (100) 

takes 

di» * „ fr, d 2 ^ ^4 a ^- (101) 
a^idi + x 2 a2 ^ic*i + ^2^2 

When the factor xi^i + 3:2^2 is multiplied into the expression in square brackets, the 
cross term is cancelled and a separated expression remains. The result is 

H a = — — [ ^(A 1 d x -{b x +x 1 d 1 ){b 2 + x 1 d 1 )x 1 d l (102) 

+A 2 d 2 - (-61 + x 2 9 2 )(-62 + x 2 d 2 )x 2 d 2 ), 



where a = b\ + b 2 , [3 = b\b 2 . This has a product of two generalized hypergeometric 
functions as a solution, and the desired result has been obtained. 

An operator of the form occuring in the three-body Toda equation can be obtained 
from ( p9|) by a few simple transformations which indicate other natural forms of this 
2-dimensional hypergeometric equation: A Fourier transform leads to an operator 
where the isolated d\, d 2 become x\ and x 2 . A point transformation X\ 1— > x\, 
x 2 1— > x\ has the mild effect X\d\ 1— > k~ l x\di, x 2 d 2 1— > k^ 1 x 2 d 2l so one can consider 
powers of the isolated factors different than 1. Finally taking k = 2 and Ai = 1 = A 2 , 
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an exponential point transformation x\ \— > e Xl } x 2 i— > e^ 2 results in an operator of the 
form that the three-body Toda equation takes after [VT2] in fl83l). 

A superficial difference in the treatment of Toda is that the Bessel functions are 
eigenfunctions with nonzero eigenvalue of the Bessel operators and not simply the 
zero-eigenvalue solutions, as the hypergeometric functions are here. This is possible 
because the factors multiplying the Bessel operators conspire to cancel the leading 
factor (xidi + a; 2 <9 2 ) _1 when they both act on the same function. Appropriate A\ 



and A 2 would allow a similar freedom in ( |102|) . Or, one could use further Gamma 



function canonical transformations to cancel A% and A 2 and substitute operators of 
appropriate form. Using this freedom can be a powerful way to obtain more explicit 
forms of solutions to the equation than found by simply using the immediate solutions 
in terms of generalized hypergeometric functions. 

As well, one can make the observation that the "eigenvalue" is implicit in the 
definition of the hypergeometric operators. By changing the constant f3 in (P5|), 
one gets solutions of the remaining operator with different eigenvalue. The Bessel 
operator itself cannot be put in hypergeometric form until the eigenvalue is specified, 
and alone it is in this sense incomplete. 

It is not difficult to further generalize ( p9|) by multiplying the term in brackets by a 
product of polynomials in x\d\, x 2 d 2 . The remaining terms must also be appropriately 
modified so that the operator which results after the transformation (|100[) has a 
separable product solution. Defining; 

B i = nc&i? ( io3 ) 



fc x =i 

S2 



B 2 = 

k 2 =l 

the appropriate generalization is 

H = B 2 A 1 d 1 + B 1 A 2 d 2 - B 1 B 2 [(x 1 d 1 ) 2 + (x 2 d 2 ) 2 -x l d l x 2 d 2 (104) 
+a(xidi - x 2 d 2 ) + 13]. 

After making the transformation induced by (|100|) , this becomes 



Ra = -[B 2 [A 1 d 1 -B 1 {b l + x 1 d 1 ){b 2 + x 1 d l )x 1 d l \ (105) 

xidi + x 2 a 2 



-5i[A 2 9 2 - B 2 (-bi + x 2 d 2 )(-b 2 + x 2 d 2 )x 2 d; 
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It is clear that the essential step in solving (|99| ) is the transformation T(x%di + 
X2<92). In the series expansion of the final solution, if the powers of x\ and x 2 are 
indexed by n\ and n 2 , respectively, this produces a factor of \jV{n\ + n 2 ). It is 
not hard to see that this kind of trick could be fruitfully applied to uncouple other 
equations and to produce coupled factors in the series expansions of other functions. 

7 Conclusion 

Quantum canonical transformations have been used to give a simple and economi- 
cal derivation of the integral representations and Kummer solutions of the confluent 
hypergeometric and hypergeometric equations. It is straightforward to extend this 
discussion to derive the differentiation formulae relating hypergeometric functions 
of different indices. The transformations to the specific forms of the equations for 
the many classical special functions are also easily handled with canonical trans- 
formations. The quantum canonical transformations are seen to provide a useful 
systematization of techniques for solving differential equations. 

Two integral representations for the physically interesting solutions of the non- 
periodic three-body Toda equation were also found using quantum canonical transfor- 
mations. These motivated a nontrivial two-dimensional generalized hypergeometric 
equation of which the non-periodic three-body Toda equation is a special case. The 
methods described here can be applied to solve other nonseparable partial differential 
equations which do not separate in the standard orthogonal coordinate systems. 
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